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Abstract 

If the low energy cross section for 3 He + p — > 4 He + e + + u e , the 'hep' reaction, 
is > 20 times larger than the best (but uncertain) theoretical estimates, then this 
reaction could significantly influence the electron energy spectrum produced by solar 
neutrino interactions and measured in the SuperKamiokande, SNO, and ICARUS 
experiments. We compare predicted energy spectra for different assumed hep fluxes 
and different neutrino oscillation scenarios with the observed SuperKamiokande 
spectrum. The spectra with enhanced hep contributions provide better fits to the 
SuperKamiokande data. 



1 Introduction 



One of the primary science goals of the SuperKamiokande [1,2], SNO [3], and 
ICARUS [4] solar neutrino experiments is to determine the shape of the so- 
lar neutrino spectrum between ~ 5 MeV and 14 MeV. In this energy range, 
neutrinos from the /3-decay of 8 B are expected, according to solar model calcu- 
lations carried out using the best available nuclear physics data, to dominate 
the solar neutrino spectrum [5]. The shape of the neutrino energy spectrum 
from a single /3-decaying source is independent of all solar physics to an accu- 
racy of 1 part in 10 5 [6]. Hence, a measurement of the shape is a direct test 
of whether something happens to the solar neutrinos after they are created, 
i. e., of the minimal standard electroweak model. 
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The SuperKamiokande collaboration has provided [1] preliminary data for the 
energy distribution of recoil electrons created by solar neutrinos scattering off 
electrons in their detector. The data are presented in 15 bins between 6.5 MeV 
and 14 MeV and one higher-energy bin, 14 to 20 MeV, for a total of 16 bins. 
The three highest energy bins show a relatively large number of events, more 
than would have been expected from the most popular neutrino oscillation 
parameters discussed prior to the first detailed report of the energy spectrum 
by the SuperKamiokande collaboration [1]. 

Could this excess of high energy events be caused by hep neutrinos, which have 
an endpoint well beyond the ~ 14 MeV endpoint of the 8 B energy spectrum? 

The hep reaction, first discussed in connection with solar neutrinos by Kuzmin [7], 

3 He + p -> 4 He + e + + u e , (1) 



produces neutrinos with an endpoint energy of 18.8 MeV, the highest energy 
expected for solar neutrinos. It was pointed out about a decade ago [8] that 
solar neutrino detectors that measure individual recoil electron energies, like 
SuperKamiokande [1], SNO [3], and ICARUS [4], might be able to detect the 
hep neutrinos. The total rate is expected to be low, but the background is 
small in this energy range. 

The goal of this paper is to demonstrate the sensitivity of inferences regarding 
the distortion of the solar neutrino energy spectrum to assumptions regarding 
the low energy cross section factor, So (hep), for the hep reaction ( Eq. [1]), 
and to emphasize the importance of experimental and theoretical studies of 
the possible contribution of hep neutrinos. 



2 Solar model calculations 



For a given solar model, the flux of hep neutrinos can be calculated accurately 
once the low-energy cross section factor for reaction (1) is specified. The rate 
of the hep reaction is so slow that it does not affect solar model calculations. 
Using the uncertainties given in Ref. [9] for the solar age, chemical composi- 
tion, luminosity, radiative opacity, diffusion rate, and for all nuclear quantities 
except So(hep), we calculate a total uncertainty in the hep flux of only 3% if 
the S— factor is known exactly. The best-estimate hep flux is very small [9]: 

0(hep) = 2.1(1.0 ± 0.03)[ S °^) H x 10 3 cnrV 1 . (2) 

^ v 11 v n 2.3 x 10- 20 keV b J v ; 
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The bracketed-factor in Eq. (2) is equal to unity for the currently-recommended 
value for So (hep) (see discussion in the following section). 

The best-estimate 8 B neutrino flux is more than 2000 times larger than the 
flux of hep neutrinos given by Eq. (2) if the bracketed factor is set equal to 
unity. This is the reason why all previous discussions of the measurement of the 
energy spectrum in the SuperKamiokande, SNO, and ICARUS experiments 
have concentrated on the recoil electrons produced by 8 B neutrinos. Even 
with the high event rate of SuperKamiokande(~ 6800 solar neutrino events 
observed in 504 days), only a few hep interactions are expected for the standard 
estimate of S (hep). 



3 Calculated hep Production Cross Sections 

Table 1 lists all the published values of the low-energy cross section factor, 
So (hep), with which we are familiar. Since reaction (1) occurs via the weak 
interaction, the cross section for hep neutrinos is too small (see Eq. 3 below) 
to be measured in the laboratory at low enough energies to be relevant to 
solar fusion and must therefore be calculated theoretically. We have also given 
in Table 1, in the column next to each cross section estimate, one or more 
characteristic features of the physics that was used to estimate the cross sec- 
tion. A review of the history of calculations of So (hep) provides insight into 
the difficulty of obtaining an accurate value. 

The first estimate of the cross section factor by Salpeter [10] considered only 
the overlap of an incoming continuum wave function with that of a bound 
nucleon in 4 He, obtaining a large value for S (hep), ~ 300 times the current 
best-estimate. However, Werntz and Brennan [11] pointed out that if one 
approximates the final 4 He state by (Is) 4 and the initial state by (ls) 3 (s c ) 
(where s c is a continuum initial state), and antisymmetrizes in space, spin, 
and isospin, then the matrix element of the usual allowed /3-decay operator 
vanishes between the initial and final states. They obtained a cross section 
factor more than two orders of magnitude smaller than the single-particle 
Salpeter estimate. 

Werntz and Brennan [11,12] derived and used in an exploratory way a sug- 
gested proportionality between the /9-decay matrix element, Mp (which can- 
not be measured), for the reaction 3 He(p,e + z/ e ) 4 He, and the neutron-capture 
matrix element (which can be measured), M 7 , for the reaction 3 He(n,7) 4 He. 
Their derivation, which was intended only to give a crude estimate of the cross 
section factor, neglected initial state interactions and the small D— state con- 
tributions of the 3 He and 4 He ground states and also assumed the dominance 
of meson exchange currents. 
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Tegner and Bargholtz [13] stressed the importance of the D-state components 
of the 3 He and 4 He wave functions and argued that the matrix elements for 
nucleon capture on 3 He are dominated by one-body operators rather than 
the two particle meson exchange terms. They derived a new proportionality 
relation between Mp and M 7 , which they used to estimate a rather large range 
of possible values for So (hep). Wolfs et al. [14] and Wervelman et al. [15] 
measured accurately the thermal neutron capture rate for 3 He(n,7) 4 He and 
used the proportionality relation of Tegner and Bargholtz to estimate values 
of S (hep). 

Carlson [16] revealed another layer of complexity by performing a detailed 
calculation with sophisticated wave functions, showing the presence in their 
model of strong destructive interference between the mesonic exchange cur- 
rents and the one-body matrix elements connecting the small components of 
the wave functions. In the most comprehensive calculation to date, Schiavilla 
et al. [17] included a more consistent treatment of the A-isobar current and 
investigated the sensitivity of So (hep) to the assumed details of the nuclear 
physics. They found a range of values 

So(hep) = (2.3 ± 0.9) x l(T 20 keV b, (3) 

which corresponds to a fusion cross section of ~ 1CT 50 cm 2 at solar ther- 
mal energies. The central value of this range, So, cent. (hep), was adopted by 
Adelberger et al. [18] and Bahcall and Pinsonneault [9] as the best available 
estimate. A value of So(hep) in the range 20 — 30 times So, CC nt.(hep) would 
be consistent (see discussion in the following section) with all the available 
evidence from solar neutrino experiments. 

Is it possible to show from first-principle physics that S (hep) cannot exceed, 
e. g., 10 times S , cen t.(hep)? We have been unable to find any such argument. 
Therefore, for the last decade we have not quoted a total uncertainty in the 
calculated standard model predictions for the hep neutrino fluxes, although 
well-defined total uncertainties are given for all of the other fluxes [8,9]. 

The reason it is difficult to place a firm upper limit is, as emphasized by Carl- 
son et al. [16] and Schiavilla et al. [17], that the calculated value of So(hep) is 
sensitive to the model used to describe both the ground sate and the continuum 
wave functions and to the detailed form of the two-body electroweak inter- 
actions. The matrix element Mp contains separate contributions from both 
the traditional single particle Gamow- Teller operator and the axial exchange- 
current operator. In the most detailed calculations [16,17], there is a delicate 
cancellation between comparable contributions from the one-body and the 
two-body operators. For example, if one artificially changes the sign of the 
principal exchange current contribution relative to the sign of the one-body 
axial current in the calculation described in Table III of Ref. [16], the size of 
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Table 1 

Calculated values of So (hep). The table lists all the 
published values with which we are familiar of the low 
energy cross section factor for the hep reaction shown 
in Eq. 1. 



bo(hep) 

Kev 0) 


Physics 


Year 


Reference 


630 


single particle 


1952 


[10] 


3.7 


forbidden; Mp oc M 7 


1967 


[11] 


8.1 


better wave function 


1973 


[12] 


4-25 


D-states + meson exchange 


1983 


[13] 


15.3 ±4.7 


measured 3 He(n, "y) 4 He 


1989 


[14] 


57 


measured 3 He(n, ^y) 4 He 


1991 


[15] 




shell model 






1.3 


destructive interference, 


1991 


[16] 




detailed wavefunctions 






1.4-3.1 


A-isobar current 


1992 


[17] 



the calculated So (hep) is increased by a factor of 32. 

For non-experts, it is instructive to compare the calculations of the pp and 
hep reactions. The pp fusion reaction [19] occurs via the allowed Gamow- 
Teller /3-decay matrix element whereas the hep transition is forbidden. For the 
pp reaction, the difficult-to-evaluate mesonic exchange corrections and matrix 
elements connecting small components of the wave functions are only small 
corrections (~ a few percent) to the total cross section. For the hep reaction, 
the exchange corrections and matrix elements involving small components of 
the wave function are the whole story. For the pp reaction, the effective range 
approximation allows one to use measured data to calculate to good accuracy 
the low energy fusion cross section. The somewhat analogous scaling laws 
relating hep fusion to the measured cross section for 3 He(n,7) 4 He reaction are 
not valid because low energy nucleon capture by 3 He occurs via competing 
and cancelling small effects and because of different initial state interactions. 
Hence, the estimated uncertainty in the low energy pp fusion cross section is 
small(~ 2% [18]), whereas the uncertainty in the hep cross section is much 
larger and is difficult to quantify. 
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4 Global fits to solar neutrino data 



We have investigated the predicted effects on solar neutrino experiments of an 
arbitrary size hep flux, which we will parameterize by multiplying S 0j ccnt (hep) 
by a constant a that is much greater than unity (cf. Eq. 3 ), 

a = ^ o(hep) (4) 

(2.3 x 10- 20 keV b) 1 ' 



We have carried out global fits to all of the solar neutrino data, the mea- 
sured total event rates in the chlorine [20], GALLEX [21], SAGE [22], and Su- 
per Kamiokande [1] experiments, the SuperKamiokande energy spectrum [1], 
and the SuperKamiokande zenith-angle dependence of the event rate [1]. We 
use the methods and the data described fully in Ref. [23], hereafter BKS98 . For 
MSW fits, the degrees of freedom (d.o.f.) are: 4 (total rates in 4 experiments) 
+15 (normalized spectrum for 16 bins) +9 (normalized angular distribution 
for 10 bins) —2 (oscillation parameters) —1 (hep flux) or 25 d.o.f. . For vacuum 
oscillations, the Day-Night asymmetry (1 d.o.f.) is a more powerful discrim- 
inant than the zenith- angle distribution [23]. Hence, for vacuum oscillations 
we have 17 d.o.f. . The only substantial difference from BKS98 is that in the 
present paper we find the best-fit to the measured SuperKamiokande energy 
spectrum by including an arbitrary amount of hep neutrinos in addition to 
the conventional 8 B flux. The contribution of the hep flux to the total event 
rates is negligible for all of the best-fit solutions. 

Figure 1 and Table 2 summarize our principal results. We see from Fig. 1 
that one can obtain good fits to the reported SuperKamiokande [1] energy 
spectrum for all three neutrino scenarios: no oscillations, MSW, and vacuum 
oscillations. 

Table 2 shows the best global fits to all the data that are possible by allowing 
large enhancements of the current best-estimate hep flux. The improvements 
are significant. 

The best-fit MSW solution improves from a confidence level (1 - P) of 7% to 
20% (for a = 26) even after accounting for the extra d.o.f. . A large range of 
values of a (< 30) give good fits. 

Figure 2 shows the allowed ranges of MSW parameters for a global solution 
with arbitrary hep flux. All three of the conventional MSW solutions [23] , small 
mixing angle(SMA), large mixing angle(LMA), and low mass(LOW) neutrino 
oscillations are allowed. In the global MSW solution with the standard hep 
flux, the LMA and LOW solutions are marginally ruled out at 99% C.L. . 
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Fig. 1. Combined 8 B plus hep energy spectrum. The total flux of hep neutrinos 
was varied to obtain the best-fit for each scenario. The figure shows the Ratio of 
the measured [1] to the calculated number of events with electron recoil energy, 
E. The measured points were reported by the SuperKamiokande collaboration at 
Neutrino 98 [1]. The calculated curves are global fits to all of the data, the chlo- 
rine [20], GALLEX [21], SAGE [22], and SuperKamiokande [1] total event rates, 
the Superkamiokande [1] energy spectrum, and the SuperKamiokande [1] Day-Night 
asymmetry. The calculations follow the precepts of BKS98 [23] for the best-fit global 
solutions for a standard 'no-oscillation' energy spectrum, as well as MSW and vac- 
uum neutrino oscillation solutions. The horizontal line at Ratio = 0.37 represents 
the ratio of the total event rate measured by SuperKamiokande to the predicted 
event rate[9] with no oscillations and only 8 B neutrinos. 

For vacuum oscillations, the value of a corresponding to the global Xmin does 
not depend strongly on Am 2 and sin 2 20 within the acceptable region. The 
improvement in the C.L. for acceptance increases from 6% to 15% when an 
arbitrary hep flux is considered. 

The best-fit global MSW solution with an arbitrary hep flux has neutrino 
parameters given by Am 2 = 5.4 x 10~ 6 eV 2 and sin 2 20 = 5.0 x 10~ 3 , which 
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Table 2 

Global fits with arbitrary hep neutrino flux. The ta- 
ble lists the best-fit enhancement parameters, a, de- 
fined by Eq. (4), for three different neutrino scenar- 
ios: no oscillations, MSW, and vacuum oscillations. 
We also list the value of Xmin f° r the global fit (25 
d.o.f. for MSW fits and 17 d.o.f. for vacuum oscilla- 
tions fits) and the confidence level P at which the 
solution is rejected, as well as the expected num- 
ber of neutrino events in the 14-16 MeV bin and 
the 16-20 MeV bin for the 504 day data set of Su- 
perKamiokande(normalized to the total number of 
observed events [1] 14 — 20 MeV). The best-fit values 
for Am 2 and sin 2 28 are given in the text. For com- 
parison, we also list the results for the best-fit global 
solutions obtained in Ref. [23] for the standard hep 
flux, i. e., a = 1.0( with one less d.o.f.). 



Neutrino a bcst Xmin p 14-16 16-20 











events 


events 


No oscillations 


26 


25.3 


0.954 


62 


14 


No oscillations 


1.0 


36.6 


0.998 


70 


6 


MSW 


25 


30.7 


0.80 


64 


12 


MSW 


1.0 


37.2 


0.93 


70 


6 


Vacuum 


30 


23.0 


0.85 


66 


10 


Vacuum 


1.0 


28.4 


0.94 


69 
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are very close to the best-fit MSW parameters [23] with the standard (much 
smaller) hep flux. For vacuum oscillations, the best-fit global solution has 
Am 2 = 7.8 x 10- n eV 2 , and sin 26 = 0.71, again similar to the neutrino 
parameters for the best-fit vacuum solution [23] with the standard hep flux. 



5 Discussion 

We have calculated global fits to all the available solar neutrino data allowing 
for an arbitrarily large hep flux. We find good fits to all the data, including the 
electron recoil energy spectrum (see Fig. 1) reported by Super Kamiokande [1]. 
The best fits are obtained for hep fluxes that are > 20 times the flux predicted 
if the best-available estimate (Eq. [3]) for the low-energy cross section factor 
for the hep reaction is used in the standard solar model calculations. We have 
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Fig. 2. Global fits: MSW solutions. The figure shows the regions in MSW parameter 
space that are consistent with the total rates observed in the four solar neutrino 
experiments (chlorine, SuperKamiokande, GALLEX, and SAGE) and with the elec- 
tron recoil energy spectrum and zenith angle distribution that are measured by 
SuperKamiokande [1]. Contours are drawn at both 95% C.L. and 99% C.L. . 

been unable to find an argument from first-principle physics that rules out 
values of the cross section factor So (hep) that are as large as required by our 
best-fit solutions that are described in Table 2. 



At first glance, these results seem discouraging. If one allows a large hep flux to 
account for the enhancement at higher energies of the measured electron recoil 
energy spectrum [1], then it would seem to be very difficult to infer anything 
fundamental about neutrino physics from the measured recoil electron energy 
spectrum. After all, to a good approximation the distortion of the spectrum 
can be represented for small distortions by a single parameter [24] and we 
are suggesting that an additional (unknown) parameter be added to the fit, 
namely, the magnitude of the hep flux. 

Fortunately, the SuperKamiokande [1,2], SNO [3], and ICARUS [4] experi- 
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merits can all test for the possible existence of a large hep flux by measuring 
the energy spectrum beyond the energy corresponding to the endpoint of the 
8 B neutrino spectrum. Table 2 shows that solutions with a large admixture 
of hep neutrinos are expected to produce appreciable numbers of events more 
energetic than 14 MeV in the 504 days of observations studied so far in the 
SuperKamiokande detector. The region beyond the endpoint energy of the 8 B 
spectrum is an excellent region in which to search for rare events since the 
background is expected to be very small between 16 to 20 MeV . 

The SNO detector should be even more sensitive than SuperKamiokande at 
the highest electron energies because the neutrino absorption cross section 
on deuterium rises more rapidly with energy than does the electron scatter- 
ing cross section and because higher energy neutrinos absorbed by deuterium 
produce higher energy electrons, whereas for v — e scattering the energy is di- 
vided almost equally between recoiling electrons and final state neutrinos [25]. 
Quantitatively, we estimate that SNO would have, depending on which neu- 
trino oscillation parameters are chosen, two to three times the rate of pro- 
duction of electrons in the 14 — 16 MeV bin if the energy resolution were the 
same in the two detectors. Moreover, the energy discrimination for SNO may 
actually be better than for SuperKamiokande, further helping in determining 
the possible contribution of hep neutrinos. 

The ratio, r, of the number of detected events in the 14 — 16 MeV bin to 
the number of detected events in the 16 — 20 MeV bin should be large if — as 
predicted by the standard solar model — the only important neutrino sources 
contributing to events in this energy region are 8 B and hep. For the best 
global fits (large a), we see from Table 2, that r satisfies for SuperKamiokande 
operating characteristics 

, . . (events : 14 - 16 MeV) 

r(gl ° bal) s (events : 16-20 MeV) * 4 > < 5) 



and for the standard So (hep), 



(events : 14-16 MeV) 
r(a=1) S (events : 16-20 MeV) * 10 ' < 6) 



Equation (5) is a prediction, valid with or without neutrino oscillations, of the 
standard solar model and can be tested with the available SuperKamiokande [1] 
data. Basically, Eq. (5) is a statement that there are no other important sources 
of high energy solar neutrino neutrinos except 8 B and hep. Equation (6) is valid 
if the current best-estimate for S (hep) is correct. 

For 504 days of SuperKamiokande operation, the best global fits predict (see 
Table 2) about 12 ± 2 events in the 16 — 20 MeV energy bin, whereas the 
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standard standard fluxes with a = 1 predict ~ 6 or 7 high energy events . 
Many more events may be required before SuperKamiokande can distinguish 
empirically between the small and large a descriptions of the energy spectrum. 

Measurements at energies below the current 6.5 MeV lower limit are also very 
important. Figure 1 shows that the best-fit vacuum solution has a small upturn 
in the spectrum at the lowest available energies. The upturn is intrinsic to this 
vacuum solution; hep neutrinos are unimportant at the lower energies. 

Solar neutrino experiments may be able to determine, after several years of op- 
eration, both the contamination (at higher energies) by hep neutrinos of the 
energy spectrum and also a strong constraint (from measurements at lower 
energies) on the allowed range of distortion parameters due to neutrino os- 
cillations. We hope that the discussion in this paper will stimulate further 
experimental and theoretical considerations of the possible effects of hep neu- 
trinos. 

We are grateful to the SuperKamiokande collaboration for making available 
at Neutrino 98 their crucially important data on the electron recoil energy 
spectrum. We are also indebted to many members of the SuperKamiokande 
collaboration for informative and stimulating conversations on the possible 
relevance of hep neutrinos for their measurements. This work is an extension 
of our earlier collaboration with A. Smirnov, to whom we are grateful for 
many insightful remarks, on BKS98. We are indebted to K. Babu, G. Beier, 
S. Bilenky, J. Carlson, D. Eisenstein, A. Gruzinov, W. Haxton, C. Kolda, H. 
Robertson, R. Shiavilla, and R. Wiring for valuable comments that improved 
an initial version of this paper. JNB and PIK acknowledge support from NSF 
grant #PHY95-13835. PIK also acknowledges support from #PHY96-05140. 
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